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Abstract. Two series of models and their yields are presented in this paper The first series 
consists of 20 Mp, models with varying initial metallicity (solar down to Z = 10^^) and rota- 
tion (Uini = — 600kms^'). The second one consists of models with an initial metallicity of 
Z = 10^^, masses between 20 and 85 Mq and average rotation velocities at these metallicities (Uin; = 
600 — 800 km s^'). The most interesting models are the models with Z = 10^^ ([Fe/H]^ —6.6). In 
the course of helium burning, carbon and oxygen are mixed into the hydrogen burning shell. This 
boosts the importance of the shell and causes a reduction of the size of the CO core. Later in the 
evolution, the hydrogen shell deepens and produces large amount of primary nitrogen. For the most 
massive models (M > 60 M©), significant mass loss occurs during the red supergiant stage. This 
mass loss is due to the surface enrichment in CNO elements via rotational and convective mixing. 

The yields of the fast rotating 20 Mq models can best reproduce (within our study) the observed 
abundances at the surface of extremely metal poor (EMP) stars. The wind of the massive models 
can reproduce the CNO abundances of the carbon-rich UMPs, in particular for the most metal poor 
star known to date, HE1327-2326. 
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INTRODUCTION 

Precise measurements of abundances of extremely metal poor (EMP) stars have recently 

been obtained by Cayrel et al. lHI], Spite et al. 0], Israelian et al. ill], These provide 

new constraints for the stellar evolution models [seeiSS]. The most striking constraint 
is the need for primary ^^N production in very low metallicity massive stars. Other 
possible constraints are an upturn of the C/0 ratio with a [C/Fe] about constant or 
slightly decreasing (with increasing metallicity) at very low metallicities, which requires 
an increase (with increasing metallicity) of oxygen yields below [Fe/H]~ -3. About one 
quarter of EMP stars are carbon rich (C-rich EMP, CEMP stars). Ryan et al. [7] propose 
a classification for these stars. They find two categories: about three quarter are main 
s-process enriched (Ba-rich) stars and one quarter are enriched with a weak component 
of s-process (Ba-normal). The two most metal poor stars known to date, HE1327-2326 
is", '9] and HE 0107-5240 [10] are both CEMP stars. These stars are believed to have 
been enriched by only one to several stars and we can therefore compare our yields 
to their observed abundances without the filter of a galactic chemical evolution model 
(GCE). In an attempt to explain the origin of the abundances observed as well as the 
metallicity trends, I computed pre-supernova evolution models of rotating single stars 
with metallicities ranging from solar metallicity down to Z = 10^ following the work 
of Meynetet al. f.11.1 . 



DESCRIPTION OF THE STELLAR MODELS 



The computer model used to calculate the stellar models is described in detail in Hirschi 
et al. lll2ll . At low metallicities the mixture of the heavy elements we adopted is the 
one used to compute the opacity tables for Weiss 95's alpha-enriched composition lfl3ll . 
The mass loss rates are described and discussed in Meynet et al. [11]. Very little was 
known about the mass loss of very low metallicity stars with a strong enrichment in 
CNO elements until recently. Vink and de Koter lll4ll study the case of WR stars but 
a crucial case, which has not been studied in detail yet, is the case of red supergiant 
stars (RSG). As we shall see later, due to rotational and convective mixing, the surface 
of the star is strongly enriched in CNO elements during the RSG stage^ Awaiting for 
future studies, it is implicitly assumed in this work (as in Meynet et al. iOb) that CNO 
elements have a significant contribution to opacities and mass loss rates. This assumption 
is supported by the possible formation of molecular lines in the RSG stage. Therefore 
the mass loss rates depend on metallicity as M ~ {Z/ZqY'-^ , where Z is the mass fraction 
of heavy elements at the surface of the star. The evolution of the models was in general 
followed until core Si-burning and the stellar yields are calculated as in Hirschi et al. 
il51 . The main characteristics of the models are presented in Table[T] More details about 
the models are presented in Hirschi [16J. 

The value of 300 kms^ used for the initial rotation velocity at solar metallicity 
corresponds to an average velocity of about 220 km s^^ on the Main Sequence (MS) 
which is very close to the average observed value [see for instance 17] . It is unfortunately 
not possible to measure the rotational velocity of very low metallicity massive stars since 
they all died a long time ago. Nevertheless, there is indirect evidence that stars with a 
lower metallicity have a higher rotation velocity. This can be due to the difficulty of 
evacuating angular momentum during the star formation, which is even more important 
at lower metallicities [see 18]. Furthermore, a very low metallicity star containing 
the same angular momentum as a solar metallicity star has a higher surface rotation 
velocity due to its smaller radius (one quarter of Z© radius for 20 Mq stars). In order 
to compare the models at different metallicities and with different initial masses, the 
ratio Uini/ Ucrit is used (see Table[l]). Ucrit is the critical velocity at which matter becomes 
gravitationally unbound. Uini/Ucrit increases only as r^^/^ for models with the same 
angular momentum (7) but lower metallicity, whereas the surface rotational velocity 
increases as r^^ (7 ~ ur). The angular momentum can be compared as well but one has 
to bear in mind that it varies significantly for models of different initial masses. Finally, 
l^ini/ 'L'crit is a good indicator for the impact of rotation on mass loss. 

In the first series of models, the aim is to scan the parameter space of rotation and 
metallicity with 20 Mq models since a 20 Mq star is not far from the average massive 
star concerning stellar yields. For this series, two initial rotational velocities were used 
at very low metallicities. The first one is the same as at solar metallicity, 300km s^^ 
The second is 500 at Z=\Q-^ ([Fe/H]~-3.6) and 600km s^^ at 2=10"^ ([Fe/H]~- 
6.6). These second values have ratios of the initial velocity to the break-up velocity, 
Uini/Dcrit, around 0.55, which is only slightly larger than the solar metallicity value 
(0.44). The 20 Mq model at Z=10"^ and with 600km s"^ has a total initial angular 
momentum 7tot = 3.3 10^^ erg s which is the same as the solar metallicity 20 Mq model 



TABLE 1. Initial parameters of the models (columns 1-5): mass, metallicity, rotation velocity [km s^ ], 
total angular momentum [10^^ erg s] and Uini/Ucrit- Total lifetime [Myr] and various masses [Mq] (7-10): 
final mass, masses of the helium and carbon-oxygen cores and the remnant mass. Total stellar yields 
(wind + SN) [Mq] for carbon (11), nitrogen (12) and oxygen (13). 
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with 300km s^^ (7tot = 3.6 10^^ erg s). So a velocity of 600 kms^^ which at first sight 
seems extremely fast, is probably the average velocity at Z=10~^. In the second series 
of models, I follow the exploratory work of Meynet et al. lllin and compute models 
at Z=10^^ with initial masses of 40, 60 and 85 Mq and initial rotational velocities of 
700, 800 and 800 km s^^ respectively. Note that, for these models as well, the initial 
total angular momentum is similar to the one contained in solar metallicity models with 
rotational velocities of 300 km s^^. Since this is the case, velocities between 600 and 
800 km s^^ are considered in this work as the average rotational velocities at these very 
low metallicities. 



EVOLUTION OF THE 20 M© MODELS 

Mass loss becomes gradually unimportant as the metallicity decreases in the 20 Mq 
models. At solar metallicity, the rotating 20 Mq model loses more than half of its mass 
and at Z = 10^^ less than 0.3% (see Table[Tl). This means that at very low metallicities, 
the dominant effect of rotation is mixing for the mass range around 20 Mq. The impact 
of rotational mixing is best pictured in the Kippenhahn diagram (see Fig. [l]). During 
hydrogen burning and the start of helium burning, mixing increases the core sizes. 
Mixing of helium above the core suppresses the intermediate convective zones linked 
to shell H-buming. So far the impact of mixing at Z = 10^^ is the same as at higher 
metallicities. However, after some time in He-burning, the mixing of primary carbon and 
oxygen into the H-burning shell is important enough to boost significantly the strength 
of the shell. As a result, the size of the helium burning core becomes and remains smaller 
than in the non-rotating model. The yield of ^^O being closely correlated with the size 
of the CO core, it is therefore reduced due to the strong mixing. At the same time carbon 
yields are increased. This produces an upturn of C/0 at very low metallicities. 
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FIGURE 1. Kippenhahn diagrams of 20Mq models at Z = 10"** with Dm = Okms"^ (left) and 
600 km s"' (right). 




FIGURE 2. Stellar yields of '^c (left), '^N (center) and "^O (right) as a function of the initial metallicity 
of the models. The solid red, dashed blue and dotted black lines represent respectively the models with 
Vini/Vc ^ 0.55 (Uini=500kms-i at Z = 10"^ and Uini=600kms-' at Z = 10"^), with Uini=300kms"i 
and without rotation. For nitrogen, the horizontal mark with C05 in the middle corresponds to the value 
deduced from the chemical evolution models of Chiappini et al. [4] . 



Stellar yields of CNO elements 

The yields of ^^Cj^N and ^^O are presented in Fig.|2land their numerical values are 
given in Table[T][see 16, for more details] . The most stringent observational constraint at 
very low Z is a very high primary ^^N production [4, .6], of the order of 0.06 Mq per star. 
In Fig.^(center), we can see that only the model at Z = 10^^ and with Uini=600km s^^ 
can reach such high values. The bulk of ^'^N is produced in the convective zone created 
by shell hydrogen burning (see Fig. 1 right). If this convective zone deepens enough to 
engulf carbon (and oxygen) rich layers, then significant amounts of primary ^^N can be 
produced (^O.OIM©). This occurs in both the non-rotating model and the fast rotating 
model but for different reasons. In the non-rotating model, it occurs due to structure 
rearrangements similar to the third dredge-up at the end of carbon burning. In the model 
with Uini =600 km s^ ^t occurs during shell helium burning because of the strong mixing 
of carbon and oxygen into the hydrogen shell burning zone. 
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FIGURE 3. Left: Kippenhahn diagrams of the 85 M© model atZ = 10"^ with Vm = SOOkms-^ Right: 
The solid lines represent the chemical composition of the wind material of the of the different models at 
Z = 10^^. The hatched areas correspond to the range of values measured at the surface of giant CEMP 
stars: HE 0107-5240, [Fe/H]~ -5.3 1 10]; CS 22949-037, [Fe/H]~ -4.0 fT9',"2^; CS 29498-043, [Fe/H]~ - 
3.5 |21]. The empty triangles |22]([Fe/H]~ —4.0) and stars ^ ([Fe/H]~ —5.4, only an upper hmit is 
given for [O/Fe]) correspond to non-evolved CEMP stars. 



TABLE 2. Initial mass (1), metallicity (2) and rotation 
velocity [kms^'] (3) and stellar wind ejected masses [Mq] 
for carbon (4), nitrogen (5) and oxygen (6). 
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Models with higher initial masses at Z = 10^ also produce large quantities of primary 
nitrogen. More computations are necessary to see over which metallicity range the large 
primary production takes place and to see whether the scatter in yields of the models 
with different masses and metallicities is compatible with the observed scatter. 

EVOLUTION OF THE MODELS AT Z = 10"^ 

Contrarily to what was initially expected from very low metallicity stars, mass loss can 
occur in massive stars ifllll . The mass loss occurs in two phases. The first phase is when 
the star reaches break-up velocities towards the end of the main sequence. Due to this 
effect stars, even metal free ones, are expected to lose about 10% of their initial masses 
for an average initial rotation. The second phase in which large mass loss can occur is 
during the RSG stage. Indeed, stars more massive than about 60 Mq at Z = 10^ become 
RSG and dredge-up CNO elements to the surface. This brings the total metallicity of 
the surface to values within an order of magnitude of solar and triggers large mass loss. 
The final masses of the models are given in Table U The case of the 85 Mq model is 



extremely interesting (see Fig.^left) since it loses more than three quarter of its initial 
mass. It even becomes a WO star. 



Wind composition and CRUMPS stars 



In Fig. |3] (right), we compare the chemical composition of the wind material with 
abundances observed in non-evolved carbon rich extremely and ultra [8] metal poor 
stars. The ejected masses of the wind material are also given in Table |2| It is very 
interesting to see that the wind material can reproduce the observed abundance in two 
ways. Either, the wind material is richer than necessary and dilution (by a factor 100 for 
example for the 40 Mq models and HE1327-2326) with the ISM is needed or the wind 
has the right enrichment (for example the 60 Mq and HEl 327-2326) and the low mass 
star could form from pure wind material. The advantage of the pure wind material is 
that it has a ratio ^^C/^^C around 5 [1 1] and it can explain Li depletion. With or without 
dilution, the wind material has the advantage that it brings the initial metallicity of the 
low mass star above the critical value for its formation ll23ill . 
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